Polar liquid crystals in two spatial dimensions: the bridge from microscopic to 

macroscopic modeling 
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Two-dimensional polar liquid crystals have been discovered recently in monolayers of anisotropic 
molecules. Here, we provide a systematic theoretical description of liquid-crystalline phases for 
polar particles in two spatial dimensions. Starting from microscopic density functional theory, 
we derive a phase-field-crystal expression for the free-energy density which involves three local 
order-parameter fields, namely the translational density, the polarization, and the nematic order 
parameter. Various coupling terms between the order-parameter fields are obtained which are in 
line with macroscopic considerations. Since the coupling constants are brought into connection 
with the molecular correlations, we establish a bridge from microscopic to macroscopic modeling. 
Our theory provides a starting point for further numerical calculations of the stability of polar 
liquid-crystalline phases and is also relevant for modeling of microswimmers which are intrinsically 
polar. 

PACS numbers: 64.70.mf, 64.70.M-, 61.30.Dk 



I. INTRODUCTION 

The study of liquid-crystalline phases formed by 
banana-shaped molecules opens the door to generate 
polar directions in a partially or completely fluid sys- 
tem due to a collective alignment of the polar axis of 
the banana-shaped (or bent-core) molecules [J. So far, 
most of the liquid-crystalline phases formed by banana- 
shaped molecules are smectic 04^] , but there have also 
been a few reports of nematic phases in this area @- 
Q. In parallel, there has been a considerable amount of 
work in Watanabe's group to generate polar nematic and 
cholesteric phases in liquid-crystalline polymers [Iol4T^ |. 
Among the polar nematic phases, a nematic phase with a 
symmetry as low as Cih (or C s ) was found [13[ confirm- 
ing earlier predictions about polar nematic phases with 
low symmetry [H| . 

About 25 years ago, there has already been an early ef- 
fort to synthesize polar nematics in systems composed of 
fairly large plate-like molecules (to avoid the flipping 
and thus to generate a lack of h — > —n symmetry, with 
h being the average preferred direction usually called the 
director [TtJ )■ About the same time, compounds com- 
posed of pyramidic molecules were synthesized with the 
same goal |l8[ , but clear-cut evidence for a polar nematic 
could not be provided in either case. This early work, 
however, triggered early mode ling in the framework of a 
Ginzburg-Landau description [19( and it was pointed out 
that phases with defects, in particular with spontaneous 
splay, should play an important role in such systems. It 
was predicted that a phase with defects would occur first 
in the vicinity of the phase transition to the polar ne- 
matic phase. 

In 2003 the group of Y. Tabe [H found a two- 
dimensional polar nematic phase in Langmuir monolayers 
using the measurements of ferroelectric response and op- 



tical investigations in a low molecular weight compound 
composed of rod-like molecules. Very recently, there were 
two additional reports on a ferroelectric response of a ne- 
matic phase in three-dimensional samples in compounds 
composed of bent-core molecules [U [22j , but it is open 
whether the ferroelectric response was due to a field- 
induced reorganization of cybotactic clusters - as sug- 
gested by the authors - or due to a bulk polar nematic 
behavior of a phase containing defects of the type out- 
lined above. 

Triggered by the reports of nematic phases in banana- 
shaped molecules, a macroscopic description of polar ne- 
mat ic p hases in three spatial dimensions was derived 
[23l |24| . It turned out that the absence of parity sym- 
metry leads in such a fluid system to a number of cross- 
coupling terms between the macroscopic polarization and 
the other hydrodynamic variables, both statically and 
in the dissipative dynamic regime. In addition, it was 
found, both for reversible as well as for irreversible dy- 
namics, that there are new cross-coupling terms not 
present in typical liquid-crystalline systems not break- 
ing parity symmetry, such as, for example, reversible dy- 
namic cross-coupling terms between flow and tempera- 
ture or concentration gradients. 

Therefore, it is of high interest to have a more mi- 
croscopic description evaluating the new cross-coupling 
terms quantitatively in order to aid synthesis of new ma- 
terials for which corresponding effects can be substan- 
tial. In this paper we start such a program using a 
phase-field-crystal (PFC) model [25H27| to analyze the 
static behavior of polar phases in two spatial dimen- 
sions. This approach can be used as a bridge from micro- 
scopic to macroscopic modeling. We will systematically 
compare the results obtained from the PFC model to 
those obtained using symmetry based approaches such 
as the Ginzburg-Landau approach, a mean-field descrip- 
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tion of phase transitions neglecting fluctuations, and the 
approach of g eneralized hydrodynamics or macroscopic 
dynamics |28(. 

While in the former only variables are taken into ac- 
count that lead to an infinite lifetime for excitations in 
the long wavelength limit, the approach of macroscopic 
dynamics also incorporates variables, which relax on a 
sufficiently long, but finite time scale in the limit of 
vanishing wave number. On realizing our program we 
strongly build on the foundations given for the static 
PFC model for nematics and other phases with orien- 
tational order in two [26| and three [27| spatial dimen- 
sions. In carrying out this program it turns out that 
it is of crucial importance for polar orientational order 
to g o beyond the Ramakrishnan-Yussouff approximation 
[29(, which is usually used in the area of PFC models. 
As a matter of fact many of the cross-coupling terms 
would not be obtained if the Ramakrishnan-Yussouff ap- 
proximation were implemented. The proposed model can 
be used as a starting point to explore phase transitions 
and interfaces for various polar liquid-crystalline sheets, 
in particular including plastic and full crystalline phases 
where the translational density shows a strong ordering. 

The paper is organized as follows: in Sec. [TTl we de- 
rive a PFC model for polar liquid crystals. Then, in 
Sec. IIII1 we discuss the relation of the two symmetry- 
based approaches with the PFC model studied in Sec. 
Hll and we show that many of the coefficients arising in 
the symmetry-based approaches can be linked to micro- 
scopic expressions via the PFC model. We finally discuss 
possible extensions of the model to more complicated sit- 
uations and give final conclusions in Sec. IIV1 



II. PHASE-FIELD-CRYSTAL MODEL FOR 
POLAR LIQUID CRYSTALS 

In general, a theory for polar liquid-crystalline phases 
can be constructed on three different levels. First of 
all, a full microscopic theory where the particle inter- 
actions and the thermodynamic conditions are the only 
input is provided by classical density functional theory 
(DFT) J35&_DFT is typically used for isotropic par- 
ticles |29l. 3^-36| but analogously holds for anisotropic 
particle interactions [361 - 140} . The second level which may 
be called mesoscopic is the phase-field approach where 
lowest-order gradients of an order-parameter field are 
considered [41|. This can be performed up to fourth- 
order gradients in order to describe a stable crystalline 
state with order-parameter oscillations leading to the 
seminal PFC model of Elder and co-workers [HQ- The 
prefactors can be brought into relation with the micro- 
scopic DFT approach [25|, EH . The PFC model has been 
extensively used to study numerically freezing and melt- 
ing phenomena on microscopic length but much larger 
(diffusive) time scales |46] - [53| . Finally, the third level are 
continuum approaches [171 . |28| . l54H56j which respect the 
basic symmetries. Here, the prefactors are phenomeno- 



logical elastic constants. PFC modeling can be used 
to assign a microscopic meaning to the prefactors thus 
linking the microscopic DFT approach to the symmetry- 
based approach. 

In this section, a PFC model for polar liquid crys- 
tals in two spatial dimensions is derived from DFT by a 
systematic gradient expansion of various coarse-grained 
order-parameter fields. As a result, we get a free-energy 
functional which involves the order-parameter fields and 
their spatial derivatives. The prefactors of various con- 
tributions are expressed as generalized moments of direct 
correlation functions in the isotropic state which provides 
a bridge between microscopic density functional theory 
and macroscopic approaches. 



A. Static free-energy functional 

We consider a two-dimensional system of N anisotropic 
particles with the center-of-mass positions r*; and orien- 
tations that are characterized by the unit vectors iii with 
i G {1, . . . , N}. To provide uniaxiality, we assume the ex- 
istence of a symmetry axis for the anisotropic particles. 
Furthermore, we assume a broken head-tail-symmetry, i. 
e., we assume polar particles. This polar system is re- 
stricted to the domain A Q M. 2 with the total area 



.4 



dr 



(1) 



and kept at a finite temperature T. The polar parti- 
cles are supposed to interact in accordance with a pre- 
scribed pair-interaction potential V(fi — rjz, u\, U2). Typ- 
ical examples include particles with an embedded dipole 
moment |57H59j modeled by a dipolar hard disk poten- 
tial, colloidal pear- like particles [gO, [6l| with correspond- 
ing excluded volume interactions, Janus particles [62ll63j 
which possess two different sides, and asymmetric brush 
polymers modeled by Gaussian segment potentials 64 1 . 
We define the one-particle density field as 



iY 



p(f,u) = / y^(f- ri)S (u - in) 



with the mean particle number density 

N 

where 



<£> = 



d N f I d N u D e -«£a=i vCfW,,*,*, 

Si 



(2) 



(3) 



(4) 



is the classical canonical average of the observable D. 
Here, we introduced the notation d n x = dxi ■ ■ • dx n for 
an arbitrary vector x and n £ N. Z denotes the classical 
canonical partition function and guarantees correct nor- 
malization such that (1) = f. Furthermore, (3 = l/(k-QT) 
is the inverse temperature with the Boltzmann constant 
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fee and S± is the unit circle. The one-particle density 
p(r,u) describes the probability density p(r,u)/p to find 
a particle with orientation u at position f. Due to the 
restriction on two spatial dimensions, the orientation 
u(ip) = (cos(p), sm(tp)) is entirely defined by the polar 
angle tp. A collective ordering of a set of particles may 
lead to a macroscopic polarization whose local direction 
can be expressed by the space-dependent dimcnsionlcss 
unit vector p(f) = u(tpa(r)), that is parametrized by a 
scalar order-parameter field <Po(r). 

Under the assumption of small anisotropics in the ori- 
entation, it is now possible to expand the one-particle 
density p{r,u) with respect to the angle ip — <po(r) be- 
tween the particular orientation u and the macroscopic 
polarization p(r) into a Fourier series. Throughout this 
paper we will assume explicitly that the preferred direc- 
tion associated with dipolar order, p, and the direction 
associated with quadrupolar order, n, are parallel. We 
will therefore use p in the following. In general, these 
two types of order can be associated with two different 
preferred directions (compare, e. g., reference [l5| ). The 
expansion with respect to orientation results in the ap- 
proximation 



p(r, + V'i(^) + P{r)(p{r) ■ «) 

+ S(f)((p(f)-uf-D) 



(5) 



where the Fourier series is truncated at second or- 
der. Here, we introduced three additional dimension- 
less order-parameter fields ipi(r), P(f), and S(r). These 
order-parameter fields arc the reduced orientationally av- 
eraged translational density 



the strength of the polarization 



(6) 



P(r) = — = / du p(r, u) (p(r) ■ u) (7) 
np J Sl 



and the nematic order parameter 



S(f) = J dup(r,u)((p(r) ■ uf - ^ 



(8) 



that measures the local degree of orientational order. The 
strength P{r) of the polarization and the director p{r) 
are modulus and orientation of the polarization P{f) = 
P(r)p(r). Note that for apolar particles |2(| P{f) = 
such that apolar particles result as a special limit from 
the present theory. 

Now we refer to microscopic density functional the- 
ory which is typically formulated for spherical systems 
[30M32( | but can also be constructed for anisotropic par- 
ticle interactions (which dates back to Onsager) |36l - l40f . 
Density functional theory establishes the existence of a 
free-energy functional J-[p(f, uj] of the one-particle den- 
sity p(f, u) which becomes minimal for the equilibrium 



density. The total functional can be split into an ideal 
rotator gas functional and an excess functional: 



T[p(f, u)] = Jid[p(r, u)] + T e * c [p{r, u)} 



(9) 



The ideal gas functional is local and nonlinear, it is ex- 
actly given by 



m<iWM = df dup(r,u){\n{K 2 p(r,u))-l) (10) 
J A JSi 

where A denotes the thermal de-Broglie- wavelength. The 
excess functional .Fexc [/?(?% u)\, on the other hand, is in 
general (i. e., for a non-vanishing V(fi — T2,u\,U2)) un- 
known and approximations are needed. However, there 
is a formally exact expression gained from a functional 
Taylor expansion in the density variations Ap(r, u) = 
p(f, u)—p around a homogeneous reference density p [30(: 

PFncWr 6)] = P^I(P) ~ E fi)] (11) 



with the n-th order contributions 
4:1 W, «)] = [drfdu C W (f, u) f[ Ap(fi , r k ) . (12) 

J A JSi i=1 

Here, c^ n '(r,u) denotes the n-particle direct correlation 
function, and the notation x = [x\, . . . , x n ) for an arbi- 
trary vector x is used. The first term on the right-hand 
side of Eq. (|11[) corresponds to n = and is an irrelevant 
constant that can be neglected. We remark that also the 
first-order term (n = 1 in Eq. (fT2"j) ) vanishes since in a ho- 
mogeneous reference state (ri , u\) must be constant 
due to translational and orientational symmetry. 

For isotropic particles, various approximations based 
on expression (jlip have been proposed. The theory of 
Ramakrishnan and Yussouff [29( keeps only second-order 
terms in the expansion. This provides a microscopic 
theory for freezing both in three [2!| and two spatial 
dimensions [65[. More refined approaches include the 
third-order term [66j with an approximate triplet direct 
correlation function (67l . |68| , but a perturbative fourth- 
order theory has never been considered. Complementary, 
non-perturbativc approaches like the recently proposed 
fundamental-measure theory for arbitrarily shaped hard 
particles (36| include direct correlation functions of arbi- 
trary order. 

We now insert the parametrization ([5]) of the one- 
particle density into Eqs. (flO)) and (|TT]) in order to ob- 
tain a free-energy functional of the order-parameter fields 
tpi(r), P{r), S(r), and p(r). First, after inserting the 
density parameterization ([S]) into the ideal gas functional 
(|10[) . we expand the logarithm and truncate the expan- 
sion of the integrand at fourth order. This order guar- 
antees stabilization of the solutions (similar to the tra- 
ditional Ginzburg-Landau theory of phase transitions). 
Performing the angular integration results in the approx- 
imation 



f3T id [p(r, u)] rj F id + 



dr /id 
A 



(13) 
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with the expansion coefficients 



(") 4S"(«)^i^f4 (17) 



with the local ideal rotator gas free-energy density 

?/> 2 P 2 ^ p2£ p4 
+ — h — 1 

2 8 4 16 

S 2 Vi^ 2 ^?5 2 P 2 S 2 S 4 
-I Z± 1- 22 1 1 

8 8 8 16 256 

and the abbreviation 

F id = 27rpA(ln(A 2 p)-l) (15) 

for a constant and therefore irrelevant term. 

Secondly, we insert the density parametrization (|5|) 
into Eq. (jllj) . We will truncate this expansion at fourth 
order. Since the n-th order direct correlation function 
c («) m Eq (jup i s no t known in general, we expand it into 
a Fourier series with respect to its orientational degrees 
of freedom. By considering the translational and rota- 
tional invariance of the direct correlation function, we 
can use the parametrization c^ n+1 ' (R, 0r, fa) with R = 
. . .,Rn), fai = ((j>Ki,---,4>Rn), and</> = (fa, . . . , (f> n ) 
for the direct correlation function c^ n+1 ^ to reduce its ori- 
entational degrees of freedoms from 2n + 2 to 2n. Here, 
the new variables are related to the previous ones by 
n - r l+ i = R l u(ip Ki ), u,i = u(ipi), faxi = <Pi- iPRi, and 
fa = ipi — <pi+i- With this parametrization, the Fourier 
expansion of the direct correlation function reads 

oo 

lj : m,j — — oo 

l<i<n I 



Next, we set A = M 2 and perform a gradient expan- 
sion (2|| H^-fz^ in the order-parameter fields. For the 
term (|12[) corresponding to n = 2, this gradient expan- 
sion is performed up to fourth order in ip 2 (r) to allow 
stable crystalline phases and up to second order in all 
other order-parameter products, where we assume that 
the highest-order gradient terms ensure stability. How- 
ever, for n = 3 and n = 4 we truncate the gradient 
expansion at first and zeroth order, respectively. This 
results in the components 



F£l[fa,P,S,p\* /df/W (18) 
Jh 2 



of the static excess free-energy functional. In this equa- 
tion, the excess free-energy densities /i™c (r) are local and 
given by 



f£l = A x i>l + A 2 (V^) 2 + ^ 3 (A^) 2 + B 1 faV • (pP) + B 2 s(p • VP - P(V -p)) 
+ B 3 (vfa ■ VS - 2(p ■ Vfa) (p • V5) - 2SVfa -((p- V)p + p(V ■ p))\ 

( 19 ) 

+ P 2 (C 1 -C a (p- Ap) -C 3 (p- V)(V-p)) +C 2 (VP) 2 + C 3 ( P - VP) 2 
+ S 2 (d 1 -±D 2 (P- Ap)) +P 2 (V5) 2 , 

/£> = PiV? + P 2 ^iP 2 + P 3 ^i5 2 + P 4 5P 2 + (PrV'i + F 2 S)P(p ■ Vfa) 

+ (2F 3 faS + F 4 P 2 + F 5 S 2 ) (p • VP) + (F 3 fa + F 6 S)P(p ■ VS) , 

/« = GiVf + G 2 ?/> 2 P 2 + G^lS 2 + G 4 faP 2 S + G 5 P 2 S 2 + G 6 P 4 + G 7 S 4 (21) 



with the coefficients 
Ai=8Mg(l), A 2 = -2M°(3), A 3 = ^M° Q (5) (22) 

in the gradient expansion in ip 2 (r), that also appear - in 
a different form - in the traditional PFC model of Elder 



and co-workers [25[. The coefficients 

S 1 =4(M?(2)-Mi 1 (2)) , (23) 

P 2 = M 2 _ 1 (2) -M\(2) , (24) 

P3 = i(M 2 _ 2 (3)+M«(3)) (25) 
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belong to the terms that contain gradients and the modu- 
lus of the polarization P(r) in first order or that describe 
the coupling between gradients in the translational den- 
sity ifii(r) and gradients in the nematic order parameter 
S(r), respectively. The following three coefficients 

G 1 =4Mj(l), 



C 2 = iMi 2 (3)-Mj(3) , 



(26) 
(27) 

G 3 = -Mi 2 (3) (28) 
appear in the gradient expansion regarding P 2 (r) and 

1 



D 1 =Mi{l), D 2 = --M^3) 



(29) 



are the coefficients of the gradient expansion in S 2 (r). 
So far, all these coefficients can also be obtained by us- 
ing the second-order Ramakrishnan-Yussouff functional 
for the excess free energy. The remaining coefficients, 
however, result from higher-order contributions in our 
functional Taylor expansion. In third order, we find for 
the homogeneous terms the coefficients 





= 32M°°, 




(30) 


E 2 


= 16(Moo U 


f 2M°j) , 


(31) 


E?, 


= 4(m 00 22 + 


2M" 2 ) , 


(32) 


Ei 


= 4(2M - 21 


+ Moo) 


(33) 



and for the terms containing a gradient we find the coef- 
ficients 



F\ = -32 (tvl^ 10 -2Mgi ! + M 



(10 

01 I ' 



F, = - 



s(m i 

+ M 



20 1 Mqj 21 - 2 M^ 2 - 2 Mjf 2 



+ MJ? + M°1 



-M°r 2 -Mjr 2 + M°i, , 



i-i 



F 4 = 16 ( M^ L ~ L — 2 MqY + M 01 



(34) 
(35) 
(36) 
(37) 
(38) 



F 5 = -4( M m " - M n -;- + Mq^ - ) , (39) 



F 6 = 2 1 2 M^ 2 " 1 - 5 - 5Mq^ + 3 M 01 



-22 



-12 



3M 2 f 2 + 2m; 



01 



(40) 
(41) 



In fourth order, we only kept homogeneous terms. The 
corresponding coefficients are 



Gi = 128M°°°, 

G 2 = 192(m 00 1 01 + M^) 

G 3 =48(m 00 202 + m00 2 ) , 



(42) 
(43) 
(44) 



G 4 = 48 ( 2 M n " o r + Mn 0n + M<& ) , (45) 



G 5 = 24 (M, 
G 6 = 48M 



G 7 - 3 M 00 2 22 



000 

111 

000 



M 



(100 



(46) 

(47) 
(48) 



All the coefficients from above are linear combinations 
of moments of the direct correlation functions. These 
moments are defined through 



Mf(a) =7 r 



2ri+l -n+ 



1 wfc**' 



~(ra+l) 



(£)■ (49) 



To shorten the notation, we introduced the abbreviations 
Mf = Mf (1) and M~* = M™ 1 ™ 2 (1, 2) and used some 
symmetry considerations that are outlined in appendix 
|A1 The moments over expansion coefficients of the direct 
correlation functions depend on the particular thermody- 
namic conditions and therefore on the parameters p and 
T. 

For stability reasons, we assume that the coefficients 
of the highest-order terms in the gradients and order- 
parameter fields are positive in the full free-energy func- 
tional. If this appears not to be the case for a certain 
system, it is necessary to take into account further terms 
of the respective order-parameter field up to the first sta- 
bilizing order. 

Eqs. (|I9[) - (j2Tj) constitute the main result of the pa- 
per: it is a systematic gradient expansion of order- 
parameter fields in the free-energy functional. The pref- 
actors are moments of various direct correlation functions 
and therefore provide the link towards microscopic cor- 
relations. This is similar in spirit to PFC models j25l l43l— 
l4rl [53j for spherical particles. 



B. Special cases of the phase-field-crystal model 

We now discuss special cases of our model. First of 
all, Eqs. (|I9p - (|2"Tj) are an extension of the excess free- 
energy density for apolar particles, that was recently pro- 
posed in reference [26j. This extension comprises a pos- 
sible polarization of liquid-crystalline particles as well as 
an enlarged functional Taylor expansion that is beyond 
the scope of the second-order (Ramakrishnan-Yussouff) 
approximation. Because of that, our free-energy func- 
tional contains a few simpler models as special cases and 
is therefore the main result of this paper. Two spe- 
cial models that are known from literature and can be 
obtained from our model by setting some of the order- 
parameter fields to zero are the traditional PFC model 
of Elder and co-workers 25| for isotropic particles with- 
out orientational degrees of freedom and the PFC model 
of Lowen [26] for apolar anisotropic liquid crystals in 
two spatial dimensions. In comparison with our free- 
energy functional, the two mentioned models base on the 
Ramakrishnan-Yussouff approximation. The traditional 
PFC model has only one order-parameter field. This is 
the translational density which corresponds to ipi{r) in 
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our model. If we set all order-parameter fields that are 
related to orientational degrees of freedom in our PFC 
model to zero, i. e., P(r) = and S(r) = 0, and neglect 
all higher-order contributions for n 2 in the functional 
Taylor expansion (jlip , then we obtain the traditional 
PFC model of Elder and co-workers. The PFC model of 
Lowcn considers anisotropic particles with one orienta- 
tional degree of freedom but no polarization. Therefore, 
this PFC model results from our model for a vanishing 
polarization P(r) = 0. Also here, we have to neglect all 
contributions (11211 for n > 2. 



III. MACROSCOPIC APPROACHES 

In this section, we investigate the bridge between 
the PFC model presented in detail in the last section 
for polar liquid crystals in two spatial dimensions and 
the symmetry-based macroscopic approaches: Ginzburg- 
Landau and generalized continuum description. In ad- 
dition, we can also compare these results obtained for 
polar liquid crystals in two spatial dimensions with those 
obtained previously for non-polar liquid crystals in two 
[26l ] as well as in three [27| spatial dimensions. 

The general PFC results of this paper have been sum- 
marized in Eqs. (|TT))) - ([2"Tj) . We first analyze the terms 
given in Eq. (|19[) . which are quadratic in the variables 
and their gradients. 

We start with terms containing only the translational 
density and its gradients in Eq. (|19|) . In the vicinity of 
the smectic-A-isotropic transition one has for the smectic 
order parameter [55| 



V>(r) = Voe +MF) 

and for the density 



(50) 



(51) 



with the average homogeneous density p (compare also 
section 6.3 of reference (5(| for a detailed discussion). 
Since the total free energy must be a good scalar, the 
smectic order parameter can enter the free energy only 
quadratically. For the lowest-order terms in the energy 
density f(f), which we define as the integrand of the free- 
energy functional 



we have [72 



F[p,P,S] = [drf 

JR 2 



-aH 2 + -&i|VVf + -& 2 |AV>| 5 



(52) 



(53) 



Comparing Eq. (f5"3"]) and the first three terms in Eq. (fTTJf , 
we can make the identifications A\, A2, and A3 with 
—u. —b\, and respectively. This situation is similar 
for non-polar nematics in three spatial dimensions [27| . 
where A\, A2, and A3 arc defined with different signs, 



and for non-polar nematics in two spatial dimensions [26J , 
where one must identify 4ir 2 pA, —4:ir 2 pB, and Air^pC 
with a, 61, and &2, respectively. 

For the terms containing only the non-polar orienta- 
tional order S(r) in Eq. (|19p. we have two contributions 
to compare to other approaches. One is spatially homo- 
geneous ~£>i and the other one is quadratic in the gra- 
dients of the orientational order ~i?2- The first contri- 
bution can be directly compared with the term "jQijQij 



in de Gennes' pioneering paper [54[. Using the structure 



Qij = S(j>iPj — \8ij) for the conventional nematic order 
parameter in two spatial dimensions, we find D\ = —A 
using the original notation of reference [54| . For the gra- 
dient terms in the Ginzburg-Landau approximation one 
has at first sight two contributions to the energy density 
just using the three-dimensional expression [54j 

L^ViQjkXViQjk) + L 2 {V i Q ik ){V :j Qjk) (54) 

for two spatial dimensions. A straightforward calculation 
shows that the two contributions are in two dimensions 
identical, however, with L\ = IL2 and thus one indepen- 
dent coefficient just as for the PFC model where one has 
the contribution ~'Z?2- 

For the terms associated exclusively with orientational 
order we have, when specialized to two spatial dimen- 
sions, in the continuum description in the energy density 

Xi(V-p) 2 4 

+ L ll ( Pi v l s) 2 



■M(V^)[fe 



K 3 (px [V xp}) 2 



(55) 



where 8j 



- 5ij — piPj is the transverse Kronecker sym- 
bol projecting onto the direction perpendicular to the 
preferred direction p(r). In Eq. (|55[) . the first line is con- 
nected to gradients of the director field p(f). It contains 
in two spatial dimensions only splay and bend and no 
twist and goes back to Frank's pioneering paper [Tt], [zll ■ 
Lines 2 and 3 in Eq. (|55|) are associated with gradients 
of the nematic modulus, S(r), and with a coupling term 
~ M between gradients of the director and gradients of 
the modulus [75|, [76| . We finally note that the gradient 
terms in Eq. (|19[) are identical to the ones given in refer- 
ence H^l, where we must identify —D21 2 in the present 
paper with 2-n 2 pE in reference [26[ . This must indeed be 
the case, since polar nematics contain the case of non- 
polar nematics as a special case in the PFC approach. 

Next, we come to the terms containing only contri- 
butions of the macroscopic polarization P(r), or equiv- 
alcntly, its magnitude (modulus) P(f) and its direction 
p(r). The term ~Ci in Eq. (|19[) is the standard quadratic 
term for a Landau expansion near, for example, the 
paraelectric- ferroelectric transition j77|. It also emerges 
when the phase transition isotropic to polar nematic is 
studied in Ginzburg-Landau approximation [X9j . The 
terms that are quadratic in gradients of P(r), i. e., the 
contributions ~ C2 and ~ C3 in Eq. (|19p . can be com- 
pared to the result of a Ginzburg-Landau approach 



(56) 
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and contain two independent contributions even in the 
isotropic phase (l9j in two spatial dimensions as is easily 
checked explicitly. 

The gradient terms for the macroscopic polarization, 
or equivalently, for its magnitude P(r) and its direction 
p(r), can also be compared to the macroscopic descrip- 
tion of polar nematics [23|, [24| . For the corresponding 
terms we have 



2 ^(V^P)(V^P) + l -K l]kl {V. lPj )(V kPl ) 



1 .,(2) 

y 



(57) 



where 6 denotes deviations from the equilibrium value, 
in particular SP = P — Pq and where the tensors are of 
the form 



1 



+ K i p l p k Sf l , 
K^piPj +K 5 5± , 



K 



(2) 
-(3) _ 



(58) 

(59) 
(60) 



Eq. (j57|) represents the analogue of the Frank orien- 
tational clastic energy (~ Kij k i) with splay and bend, 
the energy associated with gradients of the modulus (~ 

(2) 

K>j ' ) , and a cross-coupling term between gradients of the 
preferred direction to gradients of the order-parameter 
modulus (~ ^Qj-fe) ~ the analogue of the corresponding 
term in non-polar nematics [75l [78j . 

The contributions ~ C2 and ~ C3 in Eq. (|19[) are the 

PFC analo gues of the contributions ~ K^j ^ and ~ Kijki 
in Eq. (|57|) . Instead of four independent coefficients in 
the macroscopic description in two spatial dimensions, 
the PFC model gives rise to two. The contribution 
has no direct analogue in the PFC model. 

Next, we start to compare cross-coupling terms be- 
tween gradients of the variables. The discussion for the 
coupling terms between gradients of the density and gra- 
dients of the orientational order closely parallels that for 
the three-dimensional non-polar ncmatic case. In Eq. 
(flU)) . the terms of interest are proportional to B s . In 
reference [27} , these are the terms ~ _B 2 . A comparison 
of these two expressions reveals that they are identical 
in structure and that one has just to take into account 
the change in dimensionality. For spatial gradients in the 
director field coupling to spatial variations in the density 
p(r) we find in the energy density [76|, [79| 



A"(V lP )[fe + Sf kPl ] (V jPfc ) 



(61) 



By comparison with Eq. (TK)]) we find X p p = B 3 S. Finally, 
we have for the terms coupling gradients of the order- 
parameter modulus S(r) to gradients of the density [76| 



(62) 



where the second rank tensor N p is of the standard uni- 
N^piPj + NfjSjj. A comparison with 



Eq. HU) yields 2N£p = B 3 and 2N?p = -B 3 . The cou- 
pling terms listed in Eqs. (|6"Tj) and (|6"2")l exist in both 
two and three spatial dimensions. Thus, in comparison 
to the hydrodynamic description of the bulk behavior, 
which is characterized by three independent coefficients, 
we find one independent coefficient in the PFC model. 
In the framework of a Ginzburg-Landau approach using 
the orientational order parameter Qij(r) we find in the 
isotropic phase 



(63) 



and thus one independent coefficient - as has also been 
the case for the non-polar PFC model in three dimensions 
[27| as well as in two dimensions [2(|. 

The contributions ~ By and ~ B2 are containing gra- 
dients of the macroscopic polarization P(r) and couple 
to density and quadrupolar order. They are unique to 
systems with polar order, or more generally, to systems 
with broken parity symmetry, since they contain one gra- 
dient and one factor P(r). Such coupling terms are not 
possible, for example, in non-polar nematics or smectic A 
phases. The term ~B± can easily be compared with the 
macroscopic description of polar nematics given in refer- 
ence [23| | . The relevant terms from Eq. (1) of reference 
[H read 



(64) 



where Sp = p — p. We thus read off immediately that 
when comparing to the PFC model we have 2j3\p = —By 
and 2/3ip = —ByP, that is one independent coefficient 
in the PFC model and two in the macroscopic descrip- 
tion. For the term ~£?2 the situation is similar. One has 
to replace in Eq. (|64[) Sp by SS, where S(r) is the mod- 
ulus of the quadrupolar nematic order parameter with 
coefficients denoted by ^4 and ^4. Then one makes the 
identifications 2/34 = —B2 and 2^4 = B2P '. For the con- 
tribution ~£?2 we can also make easily contact with the 
Ginzburg-Landau picture. For the coupling of Pj(f) and 
its gradients to quadrupolar orientational order we obtain 
to lowest order in the Ginzburg-Landau energy density 



9%jklPl(y kQij) 



(65) 



axial form Nf- 



with g ijk i = g(SikSji + SuSjk)- This term has been given 
before for the isotropic-smectic-C* phase transition in 
liquid crystals [8(| for which the polarization Pj(r) is a 
secondary-order parameter. We note that the contribu- 
tion ~ B2 in Eq. (|19|) can be brought into a form identical 
to that of Eq. (|65[) , when it is rewritten in terms of Qij (r) 

and P{r). This shows once more the close structural 
connection between PFC modeling and the Ginzburg- 
Landau approach. 

The spatially homogeneous contributions in Eq. ([2TJ]) 
can all be interpreted in the symmetry-based framework 
as well. The term ~ Ety arises near the smectic-C*- 
isotropic phase transition [8(|: QijPiPj. The terms ^E2 
and ~ E3 can be interpreted as the density dependence 
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of the terms -~P 2 and ~ QijQij in the Landau descrip- 
tion of the polar nematic-isotropic [l9j and the non-polar 
nematic-isotropic [H3| phase transitions. Finally, the con- 
tribution ~ E\ would arise in a macroscopic description 
as a term cubic in the density variations: (Sp) 3 . Typi- 
cally, such terms are considered to be of higher order in 
a macroscopic approach. The physical interpretation of 
this term is a density dependence of the compressibility 
Most of the terms in Eq. (|20]l containing one gradient, 
namely all terms containing Fi, except for i<4, can be 
interpreted in the framework of macroscopic dynamics 
as higher-order corrections to the terms ~ /3i, ~ $4, ~ 
f3i , and ~ $4 discussed above. They correspond in this 
picture to the dependence of the coefficients Pi and Pi on 
the density changes Sp(r) and variations in the modulus 
of the quadrupolar order parameter SS(r). There is one 
exception to this picture and this is the term ~ F4 in Eq. 
(|20|l . It is also this term, which has an analogue in the 
field of the Ginzburg-Landau description of ferroelectric 
materials: 



(66) 



This nonlinear gradient term has been introduced in ref- 
erence [HI and it was demonstrated by Felix et al. [82[ 
that this term leads to qualitative changes in the phase 
diagram near the paraelcctric-fcrroclcctric transition giv- 
ing rise also to incommensurate structures. 

In Eq. (|21| . spatially homogeneous terms that are of 
fourth order in the order parameters arc presented. Most 
of them are familiar from Landau energies near phase 
transitions. The first contribution, the term ~Gi, arises 
for all isotropic-smectic phase transitions [73|, l80l . |83| 
as well as for the nematic-smectic-A and the nematic- 
smectic-C transitions fl7l. l55j : ~|-0| 4 . The contribution 
~ Gq arises near the paraelectric- ferroelectric phase tran- 
sition [73, [HJ and has also been used near the isotropic- 
polar-nematic transition [l9|: ~P 4 . The term ~ G7 is 
familiar from the non-polar nematic to isotropic [Hij and 
the smectic A to isotropic [73| transitions: ~ (QijQij) 2 ■ 
The cross-coupling term ~ G3 corresponds to an analo- 
gous term for isotropic-smectic transitions [zD, H3, l83j : 
IV" I QijQij- F° r the Ginzburg-Landau description of the 
smectic-C*-isotropic transition, the term ~ G2 arises [8(| : 
\ip\ 2 P 2 . The term ~ G5 has also an analogue at the 
smectic-C*-isotropic transition, where it has not been 
discussed before. However, for the non-polar nematic 
to isotropic phase transition in an electric field one has 
shown in reference 1841 that there are two contributions: 



XlEkE n QklQnl + X2E n E n QklQ 



ki 



(67) 



The same contributions are relevant here when the ex- 
ternal electric field is replaced by the polarization P(r). 
Finally, the term ~ G4 can be viewed as the density de- 
pendence of the term QijPiPj as it emerges near the 
isotropic-smectic-C* phase transition (80j . 



IV. CONCLUSIONS AND POSSIBLE 
EXTENSIONS 

In conclusion, we systematically derived a phase-ficld- 
crystal model for polar liquid crystals in two spatial 
dimensions from microscopic density functional theory. 
Two basic approximations are involved: first, the den- 
sity functional is approximated by a truncated functional 
Taylor expansion which we considered here up to fourth 
order. Then a generalized gradient expansion in the or- 
der parameters is performed which leads to a local free- 
energy functional. The density is parameterized by four 
order-parameter fields, namely the translational density 
■01 (f) which corresponds to the scalar phase-field variable 
in the traditional phase-ficld-crystal model, the strength 
of polarization P(r), an orientational direction given by 
a two-dimensional unit vector p, and the nematic order 
parameter S(r). In the three latter quantities, the gra- 
dient expansion is performed up to second order, while 
it is done to fourth order in ip\(f) for stability reasons. 
The traditional phase-field-crystal model [42|, |4§| and the 
recently pro pos ed phase- field-crystal model for apolar liq- 
uid crystals [26[ are recovered as special cases. The addi- 
tional terms are all in accordance with macroscopic ap- 
proaches based on symmetry considerations [28l [7a | . The 
prefactors are generalized moments of various direct cor- 
relation functions and therefore provide a bridge between 
microscopic and macroscopic approaches. 

As a general feature, we find that typically the number 
of independent coefficients for the phase-field-crystal and 
the Ginzburg-Landau approaches is the same, while in 
many cases the macroscopic hydrodynamics description 
valid inside the two-dimensional polar phase leads to a 
larger number of independent coefficients. This appears 
to be a general trend, which was also found to hold be- 
fore for the comparison of phases with three-dimensional 
non-polar orientational order [27||. In fact, it also applies 
to the two-dimensional phase-field-crystal model f or sy s- 
tems with orientational order studied in reference [26j . 

The proposed functional, as embodied in Eqs. ([19]) - 
(pH) . can be used to study phenomenologically phase 
transformations, for example, in polar nematic sheets, 
interfaces between coexisting phases [85 U87| | . and certain 
biological systems that exhibit polar order |88l. |89| . Since 
our model has more parameters, we expect even more 
complicated phase diagrams than recently numerically 
discovered in the apolar phase-field-crystal model [9(| ■ 

One could also do in principle microscopic calculations 
of the bulk phase diagram for a given intcrparticle po- 
tential V(ri — T2, u\, U2) which needs the full direct cor- 
relations of the isotropic phase as an input. The sim- 
plest idea is to neglect all direct correlation functions for 
n 3 and to rely on a second-order virial expression 91| , 
where c^(fi -r 2 ,ui,u 2 ) = e^^ 1 -^'* 1 '^) - 1, or the 
random-phase approximation for mean- field fluids [6~ij ]. 
where c (2) (r*i - f 2 ,ui,U2) = — /3V(ri - ?2,Ui,U2)- 

In a next step, the analysis can be done for Brownian 
dynamics based on dynamical density functional theory 
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[92-941, which was generalized to orientational dynam- 
ics 95| and can be used as a starting point to derive 
the order-parameter dynamics [26| . This can then be 
applied to describe the translational and orientational 
relaxation dynamics, for example, for an orientational 
glass [9o| or system exposed to a periodic driving field 
[97| . Finally, it would be interesting to generalize the 
analysis to self-pr opelled p articles which are driven along 
their orientation 98l4l00j . These particles are polar by 
definition and therefore the generalization to dynamics 
of the pr esent theo ry is mandatory to derive microscopic 
theories jlOll Il02{ for their collective swarming behav- 
ior. A dynamical theory could for example be used to 
investigate the dynamical pr opert ies of bacterial growth 
patterns of proteus mirabilis [1031 ] . 
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which implies that moments that arise from each other 
by simultaneous permutations of the elements in I, m , 
and a are equal, 



■(. . .,ai 



(A2) 



,...,7/i;,.../ n 
_ l t (,■■•? 01 j 1 ■ • ■ j a i 1 ■ ■ ■ ) i 

and the invariance of the expansion coefficients (|17[) 
against complex conjugation: 



(A3) 



The last assumption is necessary to obtain physical terms 
with real coefficients in the approximation for the ex- 
cess free-energy functional. It involves the invariance of 
(R) against simultaneous reversal of the signs of the 



-00 

-l,m 



elements in I and 



— ii,...,— l n , —mi , 
l± ,.. .,l n , mi ,. . . ,m. 



...,-m„(- R l'---' i? ") 

(Ri, . . . , R n ) , 



(A4) 



Appendix A: Symmetry considerations 

In the derivation of the approximation for the excess 
frcc-cncrgy functional, a large number of expansion co- 
efficients c["^(i?) of the direct correlation functions and 
moments M| a (a) of these expansion coefficients appear. 
To reduce their total number, we used basic symmetry 
considerations that base on four invariances of the direct 
correlation functions and showed that many of the expan- 
sion coefficients and moments are equal. This is why only 
a few moments of all possible moments for different in- 
dex combinations are present in the equations (f22"]) - (|48)) 
for the coefficients in our model. These invariances are 
the translational and rotational invariance of the direct 
correlation functions, which arc considered by an appro- 
priate parametrization c^" +1 ^(i?, <f>n, (f>) and a Fourier ex- 
pansion (fT6|) of the latter, as well as the invariance of the 
direct correlation functions concerning the renumbering 
of particles, 

^ ( • ■ • , Ti , • ■ • , Tj Ui , • ■ - , 11 j , ... ) ( Al ) 

C^ ^ ( ■ • • , Tj , ■ • ■ , Ti Uj , ■ • ■ , Ui , ... ) , 



and is equivalent to the invariance of the direct correla- 
tion functions against reflection of the system at the first 
axis of coordinates. 

When the system is apolar, the liquid-crystalline par- 
ticles have head-tail symmetry. In this case, the modulus 
P(r) of the polarization is zero and its orientation p(f) 
is not defined, while the direction n(r) associated with 
quadrupolar order still exists. Then, further symmetry 
considerations lead to the following equalities between 
expansion coefficients of the direct pair-correlation func- 
tion: 

c%{R) = ~c?l{R), 

c%{R) = cf}{R), (A5) 

c%{R) = c'il{R)- 

The consequence of these equations is, that the coeffi- 
cients B\ and B2 vanish and -B3 becomes more simple. 
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